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Abstract. We present an analysis of the subdwarf O star 
Feige 67 with a fully metal-line blanketed NLTE model at- 
mosphere based on high-resolution Far-UV (912-2000 A) and 
Near-UV (2000-3400 A) spectra from new ORFEUS II echelle 
observations and the IUE final archive. The Far-UV spectra 
are heavily blanketed by iron and nickel lines, preventing the 
detection of the stellar continuum and complicating the abun- 
dance analysis. 

Important points concerning the account for blanketing by 
millions of iron-group lines and for an accurate determination 
of iron and nickel abundances are discussed: The usage of all 
theoretically and experimentally known line opacities of the 
iron-group elements, the consideration of a broad wavelength 
range for a reliable determination of the stellar continuum flux, 
and the role of interstellar reddening. This paper outlines a 
basic approach for spectral analysis of future high-resolution 
Far-UV spectra of hot compact stars. During this study we re- 
analyzed the iron and nickel abundance of our exemplary object 
Feige 67 and confirm their extraordinarily amount. 

Key words: stars: abundances - stars: atmospheres - stars: 
evolution - stars: individual: Feige 67 - subdwarfs - Ultravio- 
let: stars 



Table 1. Observation log and properties of Feige 67. Photo- 



spheric parameters are taken from Haas (1997), except for the 
Fe and Ni abundances which have been re-analyzed in this pa- 
per 



Alias 
m v 

ORFEUS II IDs 
Observation date 
Exposure time [s] 

IUE LWR IDs (low res.) 
Observation date 
Exposure time [s] 

IUE SWP ID (high res.) 
Observation date 
Exposure time [s] 

njj j column density 
Color excess Eg_y 

Tea [K] 



BD+18 2647 
12.1 

5256.1, 5256.2 
01.12./03.12. 1996 
2954 

03400, 01570, 02944 

27.05./ 16.11. 1978, 06.01. 1979 

240, 480, 597 

20488 
20.07. 1983 
10 800 

(6 ± 1) • 10 19 /cm 2 
0.016 ± 0.005 

60 000 ±4 000 



log g [cm/s 2 ] 


5.0 




n H e/n H 


5.3 


10" 2 


nc/n H 


3.2 


10" 5 


n N /n H 


1.4 


10" 5 


no/n H 


7.3 


10" 5 


n Fe /n H 


1.4 


10" 4 


nNi/nH 


4.2 


10" 5 



1. Introduction 

With the launch of FUSE in 1999 the window to high- 
resolution Far-UV spectroscopy has been opened again. In or- 
der to prepare an adequate analysis of the upcoming data for 
hot compact stars, we studied already available high-resolution 
ORFEUS II echelle spectra using state-of-the-art NLTE (non- 
local thermodynamic equilibrium) model atmospheres with 
elaborated model atoms including the opacities of millions of 
iron and nickel lines. 

This paper pursues three aims: First we want to outline a 
basic approach for spectral analysis of future high-resolution 
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Far-UV spectra under the viewpoint of iron-group line blanket- 
ing. The second aim is a presentation of possible error sources 
and an estimation of the uncertainties of the abundances de- 
rived. The third aim is an improved determination of the iron 
and nickel abundance of our demonstration object, the sdO 
Feige 67. 

A commonly used approach to determine Fe and Ni abun- 
dances is the analysis of individual, prominent line profiles 
within a narrow wavelength range. For this purpose the ob- 
served and the synthetic spectra are normalized to the contin- 
uum. However, we will show that the detection of the stellar 
continuum can be extremely difficult or impossible because of 
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Fig. 1. The combined ORFEUS II and IUE spectra of Feige 67. The small inset figures show in more detail the overlap regions 
of the three spectra, revealing calibration problems at the blue end of the IUE-LWR spectrum. The spectra are smoothed with a 
Gaussian of FWHM = 2.0 A 



the presence of millions of lines from the iron-group elements. 
Hence, derived abundances may be afflicted with substantial 
errors. We will demonstrate that reliable analyses of Far-UV 
spectra, like from e.g. ORFEUS or FUSE, require a sufficiently 
large spectral base into the Near-UV region, where line blan- 
keting is less severe, so that the model flux can be normalized 
to the continuum in that range. For this purpose we analyzed 
the combined ORFEUS II and IUE spectra of the sdO Feige 67, 
covering the wavelength range from 900 A to 3200 A. 

Feige 67 is well suited to demonstrate the difficulties and 
possible error sources related to the Fe and Ni line blanketing 
effect, because it is known for its overabundance of the ele- 
ments (probably due to radiative levitation) and because high- 
resolution spectra are available and the stellar parameters (see 
Tab. 0) are well known ( |Haas 1997| ; |Werner et al. 1998j ). Haas 
(1997) derived the effective temperature by the ionization equi- 
libria of Fe and Ni. Our re-analysis of the Fe and Ni abundances 
results in smaller values compared to the abundances derived 
by |Haas (1997 ), however, we expect no strong back-reaction 
on the atmospheric structure which in turn might lead to a dif- 
ferent T c ff . Thus we concentrate here on a thorough analysis of 
the effect of different Fe and Ni abundances on the synthetic 
spectrum. 

In the following we describe briefly the observation and 
data reduction (Sect. 2) followed by a short outline of the model 
atmospheres used (Sect. 3). In Sect. 4 we discuss the four most 
important steps of the analysis and present the major error 
sources by means of illustrative figures. The results are sum- 
marized in Sect. 5. 



2. Observations 

The echelle spectrometer was one of the focal instruments of 
the ORFEUS telescope, which was flown on its second mis- 
sion in November/ December 1996. The wavelength range is 
900-1400 A with a spectral resolution of A/AA = 10 000. Two 
separate observations of Feige 67 were obtained with a total 
integration time of 2954 s (see Tab. [l]). The two echelle im- 
ages were co-added and then the standard extraction procedure 
was applied ( jBarnstedt et al. 1999 ). Feige 67 has been observed 
with the IUE satellite in the short and in the long wavelength 
range (see Tab. [I]). For our analysis we used the so-called "pre- 
view spectra" available in the IUE final archive. 

Any relative error in the flux calibration of the different in- 
struments with respect to each other will propagate into the er- 
ror bar of the stellar parameters to be derived, because we will 
rely on the correct overall shape of the combined spectrum. 
Fig. [I] shows the good flux calibration of all three observations 
of Feige 67, because of a good match of the adjoining spec- 
tral regions. On hand of our final model we show in Fig. g the 
strong line blanketing in the Far-UV spectral region, which re- 
sults in an overall flux depression, particularly in the degraded 
spectrum. 

3. NLTE Model atmospheres 



Based on the results of Haas (1997) we calculated a plane- 
parallel non-LTE model atmosphere with T e g — 60 000 and 
lo g g = 5.0 in r adiative and hydrostatic equilibrium. According 
to Haas (1997) the Fe and Ni abundances were set to 10 and 70 
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Fig. 2. Synthetic spectrum of the final model (including interstellar Lyman line absorption) for Feige 67, demonstrating the high 
density of Fe and Ni lines. Overplotted (white line) is the same spectrum, but degraded with a Gaussian of FWHM = 2.0 A, 
matching the smoothing of the observations shown in Fig. ^J. The dash-dotted line indicates the continuum model flux. Model 
parameters are given in Tab. Ill 
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Fig. 3. Detail of the ORFEUS II spectrum (thick line) with the final model (thin line). The flux level, which falls well below 
the continuum (dash-dotted line) is matched satisfactorily, however, the identification of individual lines is not possible because 
of inaccurate line positions in the full Kurucz line list. Vertical lines at the bottom panel mark Ni/Fe lines in the model with 
gf-values > — 2.0 (bar length proportional to log gf). The spectra are smoothed with a Gaussian of FWHM = 0.25 A 
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Fig. 4. ORFEUS II spectrum (thick line) compared to two model spectra (including reddening and Hi absorption) which are 
calculated with the full Kurucz line list in one case (thin line) and the small subset of experimentally known lines in the other 
case (thin line with circles). The full list is necessary to reproduce the spectral shape particularly below 1 100 A. All spectra are 
smoothed with a Gaussian of FWHM = 2.0 A 
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Fig. 5. The interstellar H I column density is derived from this fit to the Lya profile in the ORFEUS II spectrum. The spectra are 
smoothed with a Savitzky-Golay filter (observation: M=4, n£+n#+l=9; model: M=4, n£+n^+l=17) 



times solar, respectively. In order to calculate synthetic spectra 
with different Fe and Ni abundances, a single formal solution 
of the radiation transfer equation is performed keeping fixed 
the whole model structure. 

The computer code is based on the Accelerated Lambda 
Iteration method ( [Werner & Husfeld 1985fc |Werner 198~6| ) and 
it can handle the line blanketing of iron-group elements by a 
statistical approach using superlevels and superlines with an 



opacity sampling technique (Anderson 1985; Werner & Drei- 
zler 1999) as described in peetjen et al. (1999[ ). 

A summary of the detailed model atoms used is given in 
Tab. U Test calculations for Feige 67 have shown, that Fe V and 
Fe VI are the dominant ionization stages in those atmospheric 
layers where the iron lines are formed. 

For the calculation of opacities of iron grou p elements one 
can choose between a subset of the Kurucz (Kurucz 1991) 



line list, containing only lines known from laboratory spec- 
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Fig. 6. ORFEUS II spectrum (thick line) compared to the final model which is not attenuated by the interstellar extinction (thin 
dash-dotted line) and to final model which is attenuated by interstellar extinction for two values of E(B-V), representing the error 
margin (E(B-V)=0.011 and 0.021). The inset figures show the IUE spectra for comparison. The spectra are smoothed with a 
Gaussian of FWHM = 2.0 Aand the reddened models are normalized to the observation near 3000 A 



tra, and the full one, augmented by a vast number of theoreti- 
cally computed lines. The full line list is afflicted with uncer- 
tainties in oscillator strengths and wavelength positions, pre- 
venting the identification of individual lines when comparing 
with observed spectra (see Fig. ^|). The small list on the other 
hand causes much too less total opacity. Fig. Q demonstrates 
the high relevance of using the full line list for further analysis 
in the Far-UV range. It also demonstrates the impossibility to 
determine the stellar continuum, so that analyses with model 
atmospheres must rely on a correct flux normalization in the 
Near-UV where line blanketing is less severe. 



4. Analysis 

In the following we discuss and determine the most important 
parameters which influence the observed spectrum. Due to their 
interacting character the analysis has to be done iteratively. 



4.1. Hi column density 

The interstellar Hi column density along the line of sight is de- 
rived from the Lya profile in the ORFEUS II spectrum (Fig. |]). 
We applied different models and found no dependency of the 
derived column density on the adopted Fe and Ni abundance 
(which might be caused by strong blending of the Lya profile) 
and obtained Npjj = (6 ± 1) • 10 19 /cm 2 . The same result has 
been derived from the Lya profile in the IUE-SWP spectrum. 
Note that we did not try to normalize the model continuum to a 
putative local stellar continuum near Lya, instead we relied on 
the flux normalization in the Near-UV (see below). Therefore 
interstellar reddening has to be accounted for, otherwise one 
underestimates Njjj by approximately 0. 1 dex. This in turn (as 
reddening is calculated from Njjj) would cause an overestima- 
tion of the iron-group abundance by about a factor two. 

Another source of error can be too strong smoothing of the 
observation or the model. Even convolution with a Gaussian of 
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Fig. 7. The synthetic spectrum is normalized to the observation (thick line) near A = 3000 A to avoid the strong line blanketing in 
the Far-UV. The two thin lines denote the normalized model spectrum with the normalization factor varied by ±5%, representing 
the uncertainty in normalization. The two inset figures show details from the IUE spectra with the best fit model (thin dashed- 
dotted line) plotted additionally. All spectra are smoothed with a Gaussian of FWHM = 3.0 A 



only FWHM = 0.2 A corrupts the interstellar Lya profile be- 
cause of many blending iron group lines. Therefore we used a 
low-pass filter for noise-reduction, the Savitzky-Golay smooth- 
ing filter, as described in Press et al. (1995). 



4.2. Interstellar reddening 

Based on the determined column density, the color excess can 
be calculated to E(b-v) — 0.016 ± 0.005, using the formula 



of Groenewegen & Lamers (198S ) 



N Rl = (3.8 ±0.9) x 10 21 E {B _ 



vy 



An independent determination of the color excess could not 
be realized, because it is too small to produce a detectable 
2200 A feature in the IUE-LWR spectrum. 

In the literature the applicability of the above formula is 
debated as well as the choice of the interstellar extinction law 
itself. Hence some authors try to determine such a law adapted 
to their specific observation. Generally, however, this is impos- 
sible, so that in our case we expect a systematic error of at 
least 0. 1 dex (using the interstellar extinction law as described 



in Seaton (1979)) in addition to the statistical error, which is 
estimated to 0.1 dex from Fig. || It is clear that the neglect of 
interstellar reddening in the Far-UV range results in an under- 
estimation of the stellar flux level. This would lead to an over- 
estimation of the Fe and Ni abundances, by about a factor of 
1 .0 dex in our specific case. 



4.3. Normalization of the model to the observation 

As already mentioned above, strong line blanketing in the Far- 
UV requires a model flux normalization in the Near-UV, where 
reduced blanketing does not prevent detection of the stellar 
continuum. We determined the position of the continuum at the 
red end of the spectrum above 3000 A. In Fig. ^ the resulting 
synthetic flux has been multiplied by a factor of 0.95 and 1.05, 
respectively, representing the uncertainty of fixing the contin- 
uum in the Near-UV. A determination of the continuum only 
within the Far-UV range results in an analytical uncertainty in 
the iron and nickel abundance of about a factor two (compare 
Fig. ^ with|^). Neglecting the general slope of the spectrum in 
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Fig. 8. Detailed comparison of the final model (thin line) to the ORFEUS II (top panel) and IUE-SWP (bottom panel) spectra 
(thick lines). The other two thin lines represent models with Fe and Ni abundances varied by a factor two, representing the analysis 
error. Between 1500 A and 1600 A all models overestimate the stellar flux, perhaps due to missing opacity from other iron group 
elements which are disregarded in this study. The dashed-dotted line denotes the model continuum. All spectra smoothed with a 
Gaussian of FWHM = 3.0A 
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Fig. 9. Detail of another UV wavelength region of Feige 67, where IUE and ORFEUS II spectra overlap, showing observed spectra 
(thick line) and final model spectra (thin line). The dashed-dotted line represents the model continuum. Panel (a): ORFEUS II 
vs. final model. Panel (b): IUE-SWP vs. final model. Panel (c): IUE-SWP vs. model with Fe and Ni abundances enhanced by a 
factor two. Panel (d): IUE vs. scaled model (factor 0.93). All spectra are smoothed with a Gaussian of FWHM = 0.2 A 
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Table 2. Summary of model atoms used in our model atmo- 
sphere calculations. Numbers in brackets denote individual lev- 
els and lines used in the statistical NLTE line blanketing ap- 
proach for iron and nickel. The model atom for each chemical 
element is closed by a single level representing the highest ion- 
ization stage (not listed explicitly) 



element 


ion 


NLTE levels 


lines 




1 1 


I 


16 




JO 




He 


I 


29 




91 






1 1 


32 




1 1 s 




C 


III 


58 




231 






IV 


57 




201 




N 


III 


1 











IV 


90 




368 






V 


36 




yo 




O 


IV 


30 




99 






V 


44 




86 






VI 


52 




237 




Fc 


IV 


7 


(6472) 


25 


(1027 793) 




V 


7 


(6 179) 


25 


(793 718) 




VI 


8 


(3 137) 


33 


(340 132) 




VII 


9 


(1 195) 


39 


(86 504) 




VIII 


7 


(310) 


27 


(8 724) 


Ni 


IV 


7 


(5 514) 


25 


(949 506) 




V 


7 


(5 960) 


22 


(1006189) 




VI 


7 


(9 988) 


22 


(1 110 584) 




VII 


7 


(6686) 


18 


(688 355) 




VIII 


7 


(3 600) 


27 


(553 549) 


total 




525 


(49 041) 


1845 


(6565 054) 



the Near-UV and the Far-UV may result an additional system- 
atic error of the same magnitude. 



5. Conclusion 

In order to derive the iron and nickel abundances in the sdO star 
Feige 67 we proceeded in two major steps. At first we consid- 
ered the shape of the spectrum over a broad wavelength range. 
Here we focused on the correct consideration of all effects of 
physical relevance, i.e. large model atoms, njj j determination, 
and interstellar reddening as well as on effects of technical rel- 
evance, i.e. determination of the continuum and quality of the 
flux calibration. As a result we could demonstrate the impor- 
tant influence of these details for the second step, where the Fe 
and Ni abundances are derived from the ORFEUS II spectrum. 

We confirmed the high Fe and Ni abundances of Feige 67 
and find 3.5 and 26.5 times solar, respectively (see Tab. 
which are h o wever a factor three sma ller than the results by 
frlaas (1997| ). |Becker & Butler (1995aj |blR) derived Fe and Ni 
abundances of approximately 10 times solar, based on the IUE- 
S WP observations and using models with T c s =70 000 K. The 
reasons for this more deviating result can be i) the lower, prob- 
ably more r ealistic T e g of our model (based on the analysis of 
friaas (1997 )), ii) our use of self-consistent metal-line blanketed 
models in contrast to the Becker & Butler approach, and iii) an 
increased accuracy of our analysis using combined ORFEUS II, 
IUE-SWP and IUE-LWR observations. 

At present the accuracy of Fe and Ni abundance determina- 
tions is limited mainly by uncertainties of the atomic data avail- 
able and by the S/N ratio and resolution of the spectra. There- 
fore an error in the derived abundances of a factor of about two 
has to be accepted. For a reliable analysis of Far-UV spectra 
complementary Near-UV observations are recommended. Un- 
certainties in the knowledge of interstellar reddening enlarges 
the systematic error. 

Generally, we have demonstrated the importance of cor- 
rectly considerating iron-group line blanketing for the analysis 
of future high-resolution Far-UV spectra of hot compact stars. 
This more detailed analysis technique enables to fully benefit 
from the improved high-resolution spectra of FUSE. 



4.4. Variation of the iron and nickel abundance 

After the preparatory steps discussed above one is now in the 
position to determine the abundances of Fe an Ni by fitting 
the spectral slope. Fig. |] demonstrates the sensitivity of the 
flux distribution against the simultaneous variation of the Fe 
and Ni abundances in the model by ±0.3 dex. Detailed calcu- 
lations have shown, that modifications in the Fe abundance are 
mainly noticeable in the ranges 960-1140A andl360-1510A, 
whereas modifications in the Ni abundance occur dominantly 
in the range 1030-1330 A. The resulting abundances of the best 
fitting model are given in Tab. [l|. 

A determination of the Fe and Ni abundance by fitting the 
profile of individual lines (Fig. |9|) is rather problematic because 
of, as stated above, blending by lines with inaccurately known 
wavelength positions and oscillator strengths. 
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